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r <D ; ABSTRACT 

: 

In this article we study the galactic evolution of the LiBeB elements within 
j> ■ the framework of a detailed model of the chemical evolution of the Galaxy that 

includes galactic cosmic ray nucleosynthesis by particles accelerated in super- 
bubbles. The chemical composition of the superbubble consists of varying pro- 
portions of ISM and freshly supernova synthesized material. The observational 
trends of 6 LiBeB evolution are nicely reproduced by models in which GCR come 
from a mixture of 25% of supernova material with 75% of ISM, except for 6 Li, 
for which maybe an extra source is required at low metallicities. To account for 
q 7 Li evolution several additional sources have been considered (neutrino-induced 

. nucleosynthesis, nova outbursts, C-stars). The model fulfills the energetic re- 

c$ ■ quirements for GCR acceleration. 

> : 
X: 

_c3 ; 1. INTRODUCTION 

The evolution of the light nuclides ( 6,7 Li, 9 Be, and 10 ' n B) has been, since the early 
70's, a major topic in the studies of galactic chemical evolution. The theoretical scenario for 
the synthesis of lithium, beryllium, and boron has remained unchanged for years since the 
pioneering work of Meneguzzi, Audouze, & Reeves (1971). Recent observations (see Molaro 
et al. 1997; Boesgaard et al. 1999b, for instance) have revealed an almost linear correlation 
of both Be and B with Fe at low metallicities. This new trend seems to rule out the common 
view of galactic cosmic ray (GCR) particles being accelerated from material of the interstellar 
medium (ISM), as the latter predicts a quadratic dependence on metallicity, and it has 
motivated models of LiBeB nucleosynthesis that consider new sites for the production and 



1 Institut d'Estudis Espacials de Catalunya, Edifici Nexus, Gran Capita 2-4, 08034 Barcelona, Spain 



- 2 - 



acceleration of the GCR (Ramaty et al. 2000; Fields et al. 2000; Parizot & Drury 2000). A 
source of low-energy cosmic rays, enriched in C and O, was one of the first suggested solutions 
to the problem, mainly inspired by gamma-ray observations in Orion which, however, once 
revised appeared spurious (Bloemen et al. 1999). That has led, temporarily at least, to 
discard the low-energy component hypothesis. 

Scenarios where the composition of the GCR was identical to that of the ejecta from 
individual supernovae have also been ruled out by observations with the Cosmic Ray Isotope 
Spectrometer (CRIS) on board of the Advanced Composition Explorer (ACE) that measured 
the abundance of the radioactive isotope 59 Ni and its decay product 59 Co in current epoch 
GCR (Binns et al. 1999). 59 Ni decays by electron capture with a half life of 7.6 xlO 4 years, 
but it does not decay if it is accelerated. The abundance measurements show a very small 
amount of 59 Ni in GCR, and so there must be a minimum delay of 10 5 years between the 
explosive nucleosynthesis of the material entering in GCR and its acceleration. Such a delay, 
however, could be accounted for in models where the GCR originated in a superbubble (SB) 
region of hot, rarefied, metal-rich gas swept out by the collective effects of massive star 
winds and supernova explosions. 

In the SB scenario of Ramaty et al. (2000) LiBeB nucleosynthesis is a primary mech- 
anism (i.e., production oc O), as the composition of the GCR, that exclusively comes from 
the ejecta of supernovae inside the SB, does not change in time, and it is quite metal-rich. 
Other authors preferred a combination of a primary and a secondary (i.e., production oc 
O 2 ) mechanism. Parizot & Drury (2000) considered also a SB scenario for the accelera- 
tion of GCR, but they adopted a composition of the SB material that is a mixture of ISM 
plus a 2 or 3% of SN ejecta, according to non-magnetic SB models. On the other hand, 
Fields et al. (2000) combined two sources of GCR: one accelerated from the ISM, that has 
a secondary character, and a second component coming from SB and made of pure ejecta 
of SNe, having thus a primary behaviour. In the last two cases, the primary mechanism 
dominates at low metallicities, while the secondary one takes over when the ISM becomes 
more metal-enriched. Nevertheless, despite that it seemed that the old secondary-only GCR 
nucleosynthesis paradigm was unable to reproduce the Be and B vs Fe slope, Fields & Olive 
(1999) succeeded in doing so by using the recently observed O/Fe relationship, with a single 
slope ~ —0.35 for [O/Fe] vs. [Fe/H] at all metallicities (Israelian, Garcia Lopez, & Rebolo 
1998; Boesgaard et al. 1999a; Israelian et al. 2001). The situation about [O/Fe] evolution 
is, however, far from being clear, as stated in Fulbright & Kraft (1999). 

Most of the models for light element evolution published in recent years used analytical 
approximations for the evolution of some of the important elements in the LiBeB synthesis: 
carbon, nitrogen, oxygen, and iron (Fields & Olive 1999; Fields et al. 2000; Parizot & Drury 
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2000; Ramaty et al. 2000). Some were closed models (Vangioni-Flam et al 1998, Fields & 
Olive 1999, Fields et al. 2000), which are unable to solve the G-dwarf problem, as stated in 
Pagel (1997). 

Concerning lithium evolution, several works have recently been published, involving 
different types of possible stellar sources (Romano et al. 1999), one of which, massive AGB 
stars, seems to be recently ruled out as an important source (Ventura, D'Antona, & Mazzitelli 
2000). Casuso & Beckman (2000) have proposed a model that only considers galactic cosmic 
ray nucleosynthesis, but with a GCR rate proportional to the gas ejection rate from stars 
with masses lower than 3 M . 

In this article we analyze the evolution of the light elements by means of a detailed 
numerical model for the galactic chemical evolution (Alibes, Labay, & Canal 2001; here- 
after ALC). Our model assumes that the Galaxy builds by infall of extragalactic material, 
in an inside-out scenario for Galaxy formation, and adequately reproduces the main solar 
neighborhood observational constraints: age-metallicity relation, G-dwarf metallicity dis- 
tribution, current star formation (SFR) and supernova rates, and the evolution of the main 
elements up to the iron peak (see figures and tables in ALC). In section 2 the model is briefly 
described and the results are presented in section 3. The conclusions are given in section 4. 



2. THE MODEL 

We solve numerically the classical equations for the chemical evolution of the Galaxy, 
relaxing the instantaneous recycling approximation. The same ingredients as in ALC are 
used for our LiBeB evolution model. That includes the double infall of Chiappini, Matteucci, 
& Gratton (1997) given by: 

= A e -t/vr + B e -(t-W)/^ (1) 

dt 

where a{t) is the total surface mass density, tt and td are the characteristic infall timescales 
for the halo-thick disk phase and for the thin disk phase, respectively, and t max is the 
time of maximum accretion into the thin disk. We take t max = Tt = 1 Gyr and r D = 7 
Gyr. The coefficients A and B are fixed by imposing that they reproduce the current solar 
neighborhood values of the total and halo mass surface densities. 

An enriched infall model is adopted, which accretes primordial matter during the first 
Gyr, and then begins to incorporate slightly enriched material, with a metallicity of Z = 0.1 
Z , as suggested by the recent observations by Wakker et al. (1999) of a massive (~10 7 M ) 
metal-enriched (0.09 Z ) cloud that is now falling into the disk. We have also considered the 
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more usual assumption of a primordial composition for the whole infall process. Although 
in the first case there is a small contribution to the light element abundances from the infall, 
the final results for both types of infall are almost indistinguishable, as it will be shown later. 

The initial mass function (IMF), used is that of Kroupa, Tout, & Gilmore (1993), 
which takes different slopes for M/M < 0.5, 0.5 < M/M Q < 1.0, and M/M > 1.0. The 
range of stellar masses considered goes from 0.08 to 100 M . 

The star formation rate is from Dopita & Ryder (1994) (eq. 2) with m — 5/3, n — 1/3 
and v = 1.2: 



^(r, t) = V Xrn J^t) M0PC ^ (2) 

applied to the solar galactocentric distance of 8.5 kpc. 

To account for the CNO evolution some of the most recent and most extensive (both in 
mass and metallicity) yield calculations are used: van den Hoek & Groenewegen (1997) for 
low- and intermediate-mass stars, Woosley & Weaver (1995) 2 for type II supernovae and 
model W7 from Thielemann et al. (1996) for type la supernovae. Our results for the CNO 
evolution are shown in Figure 1 for the two infall compositions. Carbon and nitrogen follow 
the observational data, while [O/Fe] shows a gradual increase towards low metallicity with 
a slope of ~ —0.28. This value is slightly lower than the values obtained from observations 
of the O I triplet (Abia & Rebolo 1989; Israelian et al. 1998, 2001; Boesgaard et al. 1999a) 
of « -0.35. 



2.1. LiBeB Sources 



In Meneguzzi et al. (1971), only spallation reactions induced by GCR and Big Bang 
nucleosynthesis of 7 Li were considered as sources of LiBeB. Since then, several stellar mech- 
anisms have been proposed as sources of some of the light elements: the //-process (possible 
solution to the solar system boron and lithium isotopic ratios problem), nova outbursts (pos- 
sible producers of 7 Li) and C-stars (identified producers of 7 Li). All of them are included 
in our evolution code. Massive AGB stars are not included as a source of 7 Li, since recent 
works (Ventura et al. 2000) insist on their small contribution to the evolution of this isotope. 



2 It is reminded here that, as in our previous solar neighborhood model (ALC), only half the iron yields 
calculated by these authors is adopted 
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2.1.1. Big Bang Nucleosynthesis 

7 Li is the only light element that is produced in a significant amount in standard big bang 
nucleosynthesis (BBN) models. In non-standard BBN, some amount of other isotopes could 
also be produced, but the data available at low metallicities do not indicate the presence 
of any kind of abundance plateau for them. A BBN contribution of [ 7 Li] = 2.2 3 has been 
included according to Molaro (1999), which ensures the existence of the lithium or Spite 
plateau (Spite & Spite 1982) up to [Fe/H]« -1.5. Due to this BBN production of 7 Li, 
open and closed box models of Li evolution do not give the same results, since in closed box 
models all the primordial lithium is already in place when the evolution begins. Therefore, 
closed box models with the same nucleosynthetic prescriptions than open ones do produce 
an overabundance of Li at [Fe/H] = by a factor ~ 1.6. This is not the case for Be and B, 
for which open and closed box models give similar results. 

2.1.2. Galactic Cosmic Ray Nucleosynthesis 

As in Parizot & Drury (2000), we adopt a composition of the Galactic Cosmic Rays 
(A^gcr) corresponding to their being accelerated inside a superbubble, where newly synthe- 
sized material ejected by a supernova (^Q, e j) is accelerated by the shock waves generated by 
other SN and mixed with the ISM of that epoch (X^ism), the evolution of the latter being 
that shown in Figure 1. 

A^gcrW = a c] X itC] (t) + (1 - a$)XijsM.(t) (3) 

Different compositions of the GCR are used by changing the free parameter a e j in eq. 
3. For a C j = we have pure ISM CR composition and for a c - s = 1, pure SNII ejecta CR 
composition. 

The CR source spectrum considered is the so-called shock acceleration spectrum: 

q(E) oc P -^-e- E ' E ° (4) 

where p and E are the particle momentum and kinetic energy, both per nucleon, (3 = v/c 
and, for the two parameters, we have chosen s = 2.2 and Eq = 10 GeV/n. E takes into 
account the possibility that the spectrum is cut off at high energies. Since a leaky box model 
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for GCR propagation is adopted, a escape path of CR from the Galaxy of X esc = 10 g cm 2 
has been included in the propagation calculations. 

For the calculation of the production rate by spallation reactions we have used a re- 
vised version of the LiBeB code from Ramaty et al. (1997), which includes the cross 
sections of Read & Viola (1984) and takes into account all the possible reactions leading 
to any LiBeB isotope, both direct and inverse reactions (including two-step processes like 
16 0(p,x) 12 C(p,x') 7 Li) . Since GCR particles are supposed to come from the material inside a 
SB, the GCR flux, and therefore the spallation reaction production rate, is made proportional 
to the Type II supernova rate. 

We do not consider any low-energy cosmic rays (LECR), since the only observational 
evidence of this type of CR, spotted in the Orion star formation region, has been withdrawn 
(Bloemen et al. 1999). Future X-ray and gamma-ray line observations should clarify whether 
LECR actually exist. 

2.1.3. The v -process 

Woosley & Weaver (1995), in their calculations of the yields from type II supernovae, 
took into account the contribution from the neutrino-induced nucleosynthesis, i.e. the spal- 
lation reactions between the huge flux of neutrinos ejected by the explosion and the material, 
newly synthesized, in the intermediate layers of the star. This primary mechanism, which 
has not been observationally corroborated, produces mainly 7 Li and n B, so it could be quite 
important because GCR nucleosynthesis alone can not reproduce the Solar System lithium 
and boron isotopic ratios. Since these yields are quite uncertain because the flux inten- 
sity and energy spectrum of the neutrinos are not well known, a correction factor, f u , is 
introduced, its value being obtained by fitting the meteoritic boron isotopic ratio. 

2.1.4. Novae 

Nova outbursts could be important nucleosynthetic sites for those elements which have 
overproduction factors there, relative to solar, above 1000. According to the recent yields 
for CO and ONe nova outbursts from Jose & Hernanz (1998), that is the case for 7 Li (nova 
outbursts produce, during the thermonuclear runaway, 7 Be, which decays into 7 Li). We have, 
thus, adopted a 7 Li average yield per nova outburst of 1.03 xl0~ 10 M Q , taking into account 
that 30% of the nova outbursts come from ONe white dwarfs, and an outburst rate given by 
the following equation: 
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dRoutburts = D f 9 -" 

dt JM(t)+0.5 M B 

/ f(v)^(t - Tm b (1-v) ~ t coo \)dfl (5) 

where D ensures that the current nova outburst rate is ~40 yr _1 (Hatano et al. 1997), Mb 
is the binary mass, \i is the ratio of the mass of the secondary star to that of the whole 
binary system and finally f(fi) is the binary mass-ratio distribution function of Greggio & 
Renzini (1993). The time for the WD to cool enough to be able to start producing outbursts, 
t CO oi, is set to 1 Gyr. 7 Li production by nova outbursts could be confirmed by the future 
INTEGRAL mission. 



2.1.5. C-stars 



Even if supernovae (//-process) and novae are theoretically possible sources of 7 Li, 
carbon-rich stars (C-stars) with C/O > 1 give the only undeniable evidence for an stel- 
lar origin of this isotope, since they are seen to be Li-rich. The majority of them show low 
carbon isotopic ratios ( 12 C/ 13 C< 15) and are of J type (Abia & Isern 1997), with masses 
lower that 2-3 M . The mechanism by which 7 Li is produced and ejected to the ISM is 
not well established, so the only way to estimate the yield is empirical. Abia, Isern, & 
Canal (1993) statistically analyzed the contribution of C-stars to the 7 Li nucleosynthesis 
and suggested a time-dependent production rate: 



f^<S>(M)V(t-T M )dM 
f£Q(M)V(t G -TM)dM 



where M\ = 1.2 M Q and M u = 3 M . They estimated empirically a current production 
rate, from a sample of galactic C-stars, of Sf = 2 x 10~ 9 M pc _2 Gyr _1 . The authors point 
out that this 7 Li production rate can only be a lower limit and they needed to increase it 
to 6 x 10~ 8 M pc _2 Gyr _1 (Abia, Isern & Canal 1995) in order to obtain a good fit to the 
solar system data. In the present model, this rate has only to be increased up to 1.5 x 1CT 8 
MQpc^Gyr- 1 . 
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2.2. Destruction of LiBeB 

Traditionally, LiBeB isotopes are supposed to be completely destroyed in stars, so when 
they die they do not contribute to the ISM abundance. This assumption has been checked 
by means of main-sequence stellar models (Hansen & Kawaler 1994) that include all the 
nuclear reactions involved in LiBeB destruction, whose rates are taken from Caughlan & 
Fowler (1988). We have evaluated, for the whole range of stellar masses, the fraction of 
each LiBeB isotope that survives and is ejected from the star. Our results are displayed in 
Figure 2. As expected, the two lithium isotopes are the most fragile ones within the LiBeB 
group, and they are almost completely destroyed in stars over a wide range of stellar masses. 
Beryllium and both boron isotopes are more resistant, since they are destroyed at higher 
temperatures, but nevertheless the fraction of the original BeB that survives is, as in the 
case of the lithium isotopes, very small. 

These destruction terms are included in our calculations, and there are only minute dif- 
ferences in the final results as compared with those obtained assuming complete destruction. 

3. RESULTS 

We have analyzed the evolution of the light elements, taking into account all the sources 
described in section 2.1 and the restrictions that energetics puts on spectrum and composition 
of GCR. In Table 1 a summary of all the different sources considered is presented. It is 
important to notice that no other combination of the contribution from the different sources 
could reproduce all the constraints considered (Li plateau, beryllium evolution versus iron, 
and solar system boron isotopic ratio and lithium abundance). 

Table 1: SUMMARY OF THE DIFFERENT SOURCES CONSIDERED 

source isotopes produced contribution fixed by 

BBN HJ 

GCRN LiBeB 

//-process 7 Li and 11 B 

novae 7 Li 

C-stars 7 Li 



[ 7 Li]=2.2 Li plateau 

a c] = 0.25 Be vs Fe 

U = 0.29 ( n B/ lo B) 
averaged yield 

Sf = 1.5 x 10^ 8 Mepc^Gyr- 1 Li 
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3.1. Abundances Evolution 

Concerning the GCR composition, we show in the left panel of Figure 3 that for values 
of a e j (see eq. 3) near 0.25 (solid line) good fits to the observed Be and B evolution as a 
function of [Fe/H] are obtained. The x 2 test gives an interval for the possible a C j values 
of (0.12,0.41) for a level of confidence of 90%. For values of a e j near (GCR made only 
from accelerated ISM), a slope of 2 results, while, if taking a value of a C j close to 1 (GCR 
made of SNe ejecta) the contribution from GCR nucleosynthesis to LiBeB evolution at low 
metallicities is too large. The preceding means that a double origin for the composition of 
GCR must be considered, since both a pure SNII composition and a pure ISM composition 
can be rejected. Such a mixture is not well explained in non-magnetic SB models (which 
favour a C j « 0.02), unless we assume that the bulk of the GCR comes from the central part 
of the SB, where the SN ejecta material would be more dominant (Higdon et al. 1998). 
However, magnetic SB models (Tomisaka 1992) give higher values for a e j, because magnetic 
fields reduce heat conduction, and thus the evaporation of the ISM shell surrounding the SB. 
The value of a e j for these latter models is close to our best fit. In the right panel of Figure 
3, we compare the evolution, in the a e j = 0.25 case, for the two different infall compositions. 
One sees that the evolution is practically unaffected by the chemical composition of the 
accreted matter, as long as its enrichment remains moderate. 

The very recent beryllium abundance determination in the very metal-poor halo star G 
64-12 ([Fe/H]= -3.30 and [Be]= -1.1) by Primas et al. (2000) is not well fitted. As stated 
by the authors, this value would either suggest a flattening of the Be vs Fe relationship at 
low metallicities or dispersion of values in the early Galaxy. More low-metallicity data are 
needed to confirm either of these suggestions. 

Parizot (2000), also within the framework of a SB model for the acceleration of GCR, 
analyzed the proportion of pure ejecta of supernovae accelerated in superbubbles (our pa- 
rameter ct e j). He found a best value of 0.03 when adopting a GCR energy spectrum 
q(E) oc E^ 1 exp(— E/Eq) with a rather low cut-off energy of E = 500 MeV/n, that is 
supposed to mimic the results of particle acceleration in SB by a collection of weak shocks. 
Thus, our best value of a C j = 0.25 is a factor of ~ 8 larger than the result obtained by Parizot 
(2000). Such discrepancy comes in part from the higher efficiency of his weak SB spectrum. 
Had this spectrum been used instead of that in eq. (4), a larger Be production per erg 
would have been obtained, by a factor of ~ 2.2, thus reducing by the same factor our mixing 
coefficient a C y A further factor of ~ 3 is due to the fact that Parizot (2000) did not calculate 
the evolution of the actual Be/O abundance ratio, but only that of the Be/O production 
ratio by dividing the yields of Be and O per supernova. In that way, stellar astration, which 
is important for late times and high metallicities, was not taken into account. Astration 
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severely reduces the abundance of Be, but not that of O, thus lowering the Be/O abundance 
ratio at late times (Fields et al. 2001). Our model, neglecting stellar astration, produces a 
Be abundance at [Fe/H]= which is ~ 3 times larger than the value measured in the Solar 
System. Finally, the study by Parizot (2000) does not follow the detailed evolution of the 
composition of the ISM by means of a chemical evolution model. Instead, he assumes solar 
composition scaled with metallicity, i.e. O/H, which translates into overestimated C and N 
abundances, both in the ISM and in the 1 — a e j fraction of GCR accelerated from the ISM. 

Since beryllium is supposed to be produced only by GCR, and therefore its abundance 
is expected to be directly dependent on the oxygen abundance in the ISM, it is interesting 
to compare our model results with the observational data in a [Be] vs [O/H] plot (Fig. 4). 
We notice that if we use the oxygen calculated in ALC, the best value of the parameter 
a e j does not seem to be the same as that obtained from the [Be] vs [Fe/H] plot. This flaw 
was expected since our calculated oxygen evolution does not perfectly fit the Israelian et 
al. (1998, 2001) and the Boesgaard et al. (1999a) data, although our results are closer to 
them than those obtained in other recently published galactic chemical evolution models. 
However, when their relationship for [O/Fe] vs [Fe/H] is used to transform the Be-Fe plot 
to a Be-0 one, a reasonable fit is recovered for a e j = 0.25, at least for the low metallicity 
region where our results closely follow the mean trend of the data. At high metallicity we 
slightly overpredict the Boesgaard et al. (1999b) data, but the Solar System value is well 
reproduced still. 

In Figure 5, the evolution of total lithium and 6 Li are shown: in the left panel for several 
values of the parameter a C j and in the right panel for the two infall compositions and fixed a e j. 
The lithium plateau is well reproduced, as well as the upper envelope of the population I stars 
and the Solar System value. As shown in Figure 6, the increase of the lithium abundance 
for [Fe/H] > —1.5 is due to the contribution from the neutrino-induced nucleosynthesis, 
and continues by the additional Li from C-stars and nova outbursts, whose contribution 
is delayed to intermediate and high metallicities due to the long lifetimes of the low-mass 
progenitors of those objects. C-stars are the main contributors at high metallicities, as it 
is clearly seen in the Figure. The parameter a e j has little effect on the lithium evolution, 
because at low metallicities, where changes of this parameter make the GCR nucleosynthesis 
production rate change, the Li from the Big Bang and from the a + a reactions dominate. 
In the case of 6 Li, the less abundant lithium isotope and still not sufficiently measured at 
low metallicities, our model reproduces the Solar System value, but not the low-metallicity 
observations. As indicated by Ramaty et al. (2000), this may suggest the need for an extra 
source of 6 Li in the Galaxy. A possible solution to this defect of 6 Li production would be 
a higher primordial value (some non-standard BBN models produce large 6 Li abundances: 
Jedamzik & Rehm 2001), but that would mean the existence of plateaus at low metallicities 
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also for Be and B (or a 7 Li primordial value not compatible with the Spite plateau), which 
have not been observed up to now. 

As it happens for Be and B, the Li evolution is almost independent of the assumed 
composition of the infalling matter. 

3.2. Ratios Evolution 

Isotopic and elemental ratios are crucial to analyze the importance of the different 
contributions to LiBeB abundances. The boron isotopic ratio shows, for a given value of a e j, 
the importance of the neutrino induced nucleosynthesis. This source is necessary to reach 
the solar system value. On the other hand, the boron over beryllium ratio data are quite 
constant at all metallicities, implying that the //-process can not be dominant. The data 
on the next element ratio considered, Li/B, show a steep descent from the large value at 
low metallicities, due to the Li plateau, down to the solar value. Finally, there are quite 
few available data for 6 Li, so the lithium isotopic data are not a strong constraint, the solar 
system data excepted. 

In all cases we have used a value of the parameter that restricts the contribution of the 
neutrino induced nucleosynthesis of f v = 0.29 (similar to that found by Ramaty et al. 2000), 
which allows to reproduce the Solar System boron isotopic ratio, as displayed in the upper 
panel of Figure 7. Higher values of /„ would increase the ratio 11 B/ 10 B up to values not 
allowed by the Solar System data. We see that the ratio is quite constant during the halo 
phase (a higher ratio for lower values of a e j) due to the contribution of the //-process and the 
primary character of GCR nucleosynthesis in those early epochs, and it decreases towards the 
Solar System value when the secondary behaviour of the GCR becomes important enough. 
The lack of data at low metallicity hinders a more firm conclusion, but the descending trend 
after [Fe/H] = is compatible with the current ISM ratio (3.4±0.7, Lambert et al. 1998). 

A similar behavior is shown by the B/Be ratio (lower panel of Fig. 7). The constant 
ratio is higher for lower values of a e j, and the standard GCR nucleosynthesis also dominates 
at high metallicity. 

Finally, the two other ratios considered include lithium (Figure 8, with 6 Li/ 7 Li in the 
upper panel and Li/B in the lower one). As we have already pointed out, our model can 
not account for the 6 Li at low metallicities, so the lithium isotopic ratio obtained is much 
lower than the nonzero 6 Li/ 7 Li data around [Fe/H] = —2.5. The values at intermediate 
metallicities are not reproduced, since the stars where they have been measured may have 
suffered some amount of lithium destruction (more important in the case of 6 Li than for 7 Li). 



For the Li/B ratio this lack of 6 Li is not important and our models fit well the available data. 



3.3. Energetics 

Since supernovae, either individually or in SB, are supposed to be the site of GCR 
acceleration, the total energy of a SN explosion must be larger than the energy of the GCR 
particles produced by it. The mean value of the Be/Fe ratio in low-metallicity stars measured 
by Molaro et al. (1997) and Boesgaard et al. (1999b) is about 1.37 x 1CT 6 and the average 
yield of 56 Fe from SN, according to ALC prescriptions, is QsN(Fe) = 0.0515 M = 1.1 x 10 54 
atoms. Then, as shown by Ramaty et al. (1997), one can evaluate the number of atoms of 
Be produced by each SN event: 

Be 

Q SN (Be) = — Q SN (Fe) ~ 1.5 x 10 48 atoms (7) 
re 

We have also evaluated the number of atoms of beryllium produced per erg (Q(Be)/W, 
upper panel in Figure 9), so we can show in the lower panel the amount of energy per SN 
required to accelerate GCR, as 

^sn = ^SN(Be) 
GCR Q(Be)/W 

Since the early Galaxy was poor on metals, more energy per SN was required to account 
for the Be/Fe ratio than in the present day, when the secondary mechanism also contribute 
to the synthesis of Be. In the halo phase, the value of VFgcr is about 2xl0 50 erg, so the 
energy constraint is fulfilled because only 15-20% of the kinetic energy of any supernova has 
to be involved in GCR acceleration. 



4. CONCLUSIONS 

The results of a recent and successful galactic chemical evolution model are applied to 
the particular problem of the light elements. To calculate the LiBeB evolution, it is necessary 
to know which are the abundances of hydrogen, helium, and the CNO nuclei in the ISM, 
since they are the main targets of the GCR to produce lithium, beryllium, and boron. The 
evolution of iron is also important, because traditionally the light element evolution has 
been plotted in [LiBeB] vs [Fe/H] graphs and iron has been taken as representative of the 
metallicity of the stars where LiBeB have been measured. In order to achieve this goal, we 
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have used some of the latest and most complete calculations of the yields from different types 
of stellar scenarios for nucleosynthesis: low- and intermediate-mass stars, Type la and II 
supernovae, and nova outbursts. 

Our calculations show that a mixture of 25% of SN ejecta (with upper and lower limits 
of 12% and 41%, respectively) and 75% of ISM, for the GCR composition, can account for 
the linear trend of BeB vs Fe. If the O/Fe vs Fe relationship of Boesgaard et al. (1999a) 
(somewhat steeper than the one obtained by ALC) is used, an acceptable fit to the Be vs O 
plot is also obtained. 

On the other side, the lithium plateau is easily reproduced with a BBN contribution, as 
usual, and the increase of the Li abundance after the halo phase is quite nicely reproduced 
as well by taking into account several sources (the //-process, nova outbursts, and C-stars). 

Concerning the isotopic evolution, our model fails to reproduce the early-time 6 Li/ 7 Li 
ratio, due to a lack of 6 Li at low metallicities. We point out the necessity of an extra source 
for this isotope. The other calculated ratios, however, nicely reproduce the data. 

Finally, the energy problem has been addressed, in order to check that not too large a 
fraction of the supernova energy is required to produce the light elements. 

This work has been supported by the DGESIC grant PB98-1183-C03-01 and the DGI 
grant AYA2000-0983. We thank late Dr. R. Ramaty for kindly providing us with his 
numerical code for cosmic-ray induced nucleosynthesis of the light elements. 
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Fig. 1. — Carbon, nitrogen, and oxygen evolution, from ALC. Solid line: enriched model; 
dashed line: primordial model. 
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Fig. 2. — Astration model results. Mass fraction of a LiBeB isotope that has survived and 
has been ejected from each star (/j(M)). 
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Fig. 3. — Beryllium and boron evolution. Left panel: evolution for the enriched model and 
different values of a c - y Right panel: evolution for the a e j = 0.25 model and the two infall 
compositions. 
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[O/H] 

Fig. 4. — Beryllium vs oxygen evolution for the a C j = 0.25 case. Oxygen resulting from our 
model is used in thin line plots, while oxygen from Boesgaard et al. (1999a) is used in the 
thick line ones. Data: Boesgaard et al. (1999b). 
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Fig. 5. — Lithium and Li evolution. Left panel: enriched model and same values of a C j as 
in figure 3. Right panel: both infall models for a C j = 0.25. 
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Fig. 6. — Different contributions for the lithium evolution, for the enriched infall composition 
and a c \ = 0.25. 
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Fig. 7. — 11 B/ 10 B and B/Be ratios evolution. Same values of a C j as in figure 3, for the 
enriched infall composition. 
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Fig. 8. — 6 Li/ 7 Li and Li/B ratio evolution. Same values of a C j as in figure 3, for the enriched 
infall composition. 
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Fig. 9. — Upper panel: Q(Be)/W vs [Fe/H] inoura C j = 0.25 model. Lower panel: Energy per 
supernova needed to accelerate GCR, according to our evaluation of the atoms of beryllium 
produced by GCR for each SN event. 



